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SUMMARY 


A simplified analysis has been made of the factors affecting the 
loss in lift and the shift in aerodynamic center of a swept wing due 
to its distortion under aerodynamic load. The manner in which these 
particular aeroelastic effects influence the longitudinal stability of 
an airplane has been considered. 

The results show that large variations in the aeroelastic effects 
associated with wing bending are produced by changes in aspect ratio , 
sweep angle, and thickness ratio. These variations are, in general, 
larger than those produced in wing-bending stress and with the thick- 
ness ratios being contemplated today restrict the aspect ratio to low 
values for large angles of sweep and, conversely, restrict the sweep 
angles to low values for large aspect ratios. Expressions obtained 
for the ratio of angle-of-attack change due to wing torsion to angle- 
of -attack change due to wing bending show that torsion effects tend to 
alleviate the effects of bending, but the angle-of-attack changes due 
to torsion are much smaller than those due to bending except for wings 
with low values of sweep, aspect ratio, or, in particular, a combina- 
tion of the two. Decreasing the plan-form taper ratio is another 
means for extending the combinations of sweep and aspect ratio without 
an increase in wing structural weight. When aeroelastic considerations 
are important, a more rapid increase in the structural weight of wings 
with increase in airplane size appears to occur than for designs based 
on stress alone. The choice of steel or duralumin as a structural 
material is not significant insofar as aeroelastic considerations are 
concerned. 

Some alleviation of aeroelastic effects occurs in maneuvers because 
of the inertia of the wing. In addition, the effect of the wing- 
aerodynamic-center shift on the ' longitudinal stability of an airplane 
with a horizontal tail may be alleviated to some extent by the com- 
pensating effect of the wing loss in lift provided the percentage 
reduction in wing- lift- curve slope is greater than that of the tail- 
lift-curve slope. 
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INTRODUCTION 


One of the primary objectives in airplane design is the attainment 
of efficient high-speed performance. This objective, however, must be 
compromised in order to meet other requirements such as satisfactory 
stability and control characteristics and structural integrity. In 
recent years the trend has been toward the use of swept wings as a means 
for obtaining efficient performance at transonic speeds, but the use of 
swept wings has introduced a structural design problem not inherent in 
unswept wings. This structural design problem is a consequence of the 
necessity for restricting the changes in longitudinal stability of a 
swept-wing airplane caused by the wing bending under load. 

The importance of this aeroelastic phenomenon has been recognized 
for some time, and several investigators have developed more or less 
refined methods for analysis of wing elastic distortions and their 
effects. (For examples, see refs. 1 to It.) More specific evidence as 
to the influence of wing external geometry (in conjunction with flight 
condition) on this flexibility — longitudinal-stability problem would 
appear to be of interest to the aerodynamicist, however, so that he 
can obtain rough indications as to which wing plan forms would be 
practical from aeroelastic considerations. Accordingly, a simplified 
analysis has been made of the factors affecting the loss in lift and 
the shift in aerodynamic center of a swept wing caused by itB distor- 
tion under load. A shell type of wing structure is assumed and the 
method of analysis obtains direct algebraic solutions for the extreme- 
fiber stress a n d the aeroelastic effects in terms of the wing structural 
weight and material, its external geometry, and the flight condition 
under which it operates. Initially, the effects of wing bending alone 
are considered and the effects of torsion are included subsequently 
through determination of the ratio of angle-of -attack changes resulting 
from wing torsion to those resulting from wing bending. The solutions 
are presented in the form of parametric charts which, in turn, are used 
to calculate specific examples illustrating effects produced by such 
factors as aspect ratio, thickness ratio, taper ratio, sweep angle, 
airspeed, altitude, material density, and stiffness. 

Because of the simplifications involved in the analysis, the results 
are to be considered chiefly qualitative in nature and are intended to 
be of use primarily in selecting practical wings for research programs 
or in other instances where the weighing of aeroelastic effects In the 
determination of practical ranges for wing parameters is desirable. 
Quantitative evaluation of the effects of changes in wing geometry, 
however, are believed to be provided by the results and absolute values 
of the aeroelastic characteristics of a wing may be obtained if the 
aeroelastic Characteristics of a related practical wing are known and 
this wing is used as a basis for comparison. 



SYMBOIS 


area of solid airfoil section measured perpendicular to 
elastic axis, sq ft 

. nondime ns ional area of solid airfoil section measured 
perpendicular to elastic axis, a s /T'(c') 2 

cross-sectional area of skin of assumed hollow airfoil sec 

tion measured perpendicular to elastic axis, 2 a 0 T , r(c') 
sq ft 

area enclosed by mean line of skin of airfoil section 
measured perpendicular to elastic axis, sq ft 

aspect ratio, (2s) 2 /s 

chord of swept wing in direction of free stream, ft 

chord of transformed wing perpendicular to elastic 
axis, c cos 7, ft 

mean aerodynamic chord of swept wing in direction of free 
stream, ft 

f 

wing lift coefficient, L/qS 

incremental lift coefficient due to wing bending, AT, / qS 

incremental pitching-moment coefficient about aerodynamic 
center of rigid wing due to wing bending, /M/qSc 

streamwise distance from aerodynamic center of a given 
spanwise section to aerodynamic center of rigid wing, ft 

distance from center of pressure of airfoil section to 
elastic axis, percent of chord 

Young's modulus of elasticity of wing structural material, 
lb/sq ft 

gravitational acceleration, ft/sec 2 

torsional modulus of rigidity of wing structural material, 
lb/sq ft 
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bending moment of inertia of solid airfoil section measured 
perpendicular to elastic axis, ft^ 

nondimensional moment of inertia of solid airfoil section- 
measured perpendicular to elastic axis, I S /(T' )^(C 

moment of inertia of assumed hollow airfoil section measured 

perpendicular to elastic axis, I 0 (T')^(c' ) 4 E - (1 - 2r)3, 
ft h 

torsional moment of inertia of hollow airfoil section 
measured perpendicular to elastic axis, ft^ 

chordwise location on transformed wing of front shear web 
from leading edge and of rear shear web from trailing 
edge, fraction of chord 

proportionality constant related to spanwise distribution 
of loading over wing 

tail length, streamwise distance between point for neutral 
static stability with respect to angle of attack and 
center of pressure of horizontal tail, ft 

wing lift, lb 

incremental lift due to win g bending, lb 
moment, ft-lb 

incremental, pitching moment, produced about aerodynamic 
center of rigid wing, due to wing bending, ft-lb 

normal acceleration, g units 

magnitude of uniformly distributed load, lb/sq ft 
dynamic pressure, pV^/ 2, lb/sq ft 

ratio of wing-skin thickness to airfoil-section thickness 
(assumed constant along chord) 

semispan of swept wing, ft 

semispan of transformed wing, s/cos 7, ft 
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tail area, sq ft 

t airfoil skin thickness, ft 

T ratio of airfoil-section maximum thickness to chord 

measured in direction of free stream 

T* ratio of airfoil-section maximum thickness to chord 

measured perpendicular to elastic axis, T/cos 7 

V airspeed, ft/sec 

w weight density of wing structural material, lh/cu ft 

Wq airplane gross weight, lb 

W wing structural weight, lb 

W w wing total weight, lb 

x' chordwise station on transformed wing measured from 

leading edge, positive forward, ft 

x streamwise location of aerodynamic center measured from 

leading edge of mean aerodynamic chord, positive 

forward, ft 

hx , shift in wing aerodynamic-center location due to wing 

bending, positive forward, ft 

x n streamwise location of center-of-gravity position for 

neutral static stability with respect to angle of attack 
measured from leading edge of mean aerodynamic chord, 
positive forward, ft 

( X cp)aL center of pressure of loss in lift due to wing bending, 
measured from leading edge of mean aerodynamic chord, 
positive forward, ft 

y spanwise station on swept wing (see fig. l), ft 

y' spanwise station on transformed wing (see fig. l) 

y/cos 7 , ft 

y • spanwise location of mean aerodynamic chord of swept 

wing (see fig. l), ft 
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a 


angle of attack, radians 


Ax change in local angle of attack due to aeroelastic 

distortion, radians 

7 angle of sweep of elastic axis, deg 

6 local torsional deflection of chord of transformed wing 

relative to root chord, radians 

X plan-form taper ratio, ratio of tip chord to root chord 

A angJe of sweep of axis which is assumed to contain aerodynamic 

centers of all sections along span, deg 

I floating spanwise coordinate on transformed wing 

(see fig. l), ft 

p air density, slugs/cu ft 

a wing stress in extreme fiber, lb/sq ft 

t taper ratio of thickness ratio, T tip| T root 

dz’/cLy' slope of elastic axis measured in plane perpendicular 
to chord plane of wing 


de/da 


( C V) B 


f ;j_ ( X,t ) 


f 2 (Vr) 


f 3 (X) 


variation of average angle of downwash at tail with angle 
of attack, per radian 

variation of lift coefficient with angle of attack, 
per radian 

variation of lift coefficient with angle of attack 
associated with angle-of-attack changes due to wing 
bending, per radian 

function of taper ratio related to incremental lift 
parameter 

function of taper ratio related to incremental moment parameter 

function of taper ratio defined by 1 + — + 7 — 

3 (1 + X ) 2 
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Subscripts: 

y spanwise station y 

y' spanwise station y' 

^ spanwise coordinate i 

r wing root 

wb wing-body combination 

w wing 

t tail 

B bending 

T torsion 

R rigid wing 

F flexible wing 

i inertia 

a aerodynamic 

av average 

1 first extreme loading condition (see appendix) 

2 second extreme loading condition (see appendix) 

ASSUMPTIONS AND LIMITATIONS OF ANALYSIS 


Because application of refined methods of aeroelastic analysis 
would be extremely laborious when applied to the determination of the 
effects of wide variations in wing external geometry, the method out- 
lined herein for analysis of the effect of wing distortion on longi- 
tudinal stability involves simplifying assumptions and l imi tations - 
Although these simplifications prevent the attainment of precise quan- 
titative conclusions, the objective of qualitatively indicating prac- 
tical ranges for various design parameters is believed to be attained. 
Further justification for the simplifications derives from the fact 
that many of the parameters necessary for a more refined analysis have 
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not been and possibly cannot be accurately determined in the transonic 
speed range. 

The present analysis applies to shell-wing structures (that is, 
the load is taken primarily in the skin) in which the chordwise varia- 
tion of skin t hickn ess is proportional to the local airfoil-section 
thickness, anfl the spanwise variation in skin thickness is of a char- 
acter to give a constant spanwise stress in the extreme fiber under a 
uniformly distributed loading. This structure was chosen for ease in 
computation and, although not a typical practical, case, represents an 
efficient structure because the structural material is concentrated far 
from the elastic axis. The airfoil shape is arbitrary. 

Initially, in the present analysis, the effects of wing bending 
alone are considered and, in general, the results which are presented 
were obtained for this case. The effects of torsion are included by 
subsequent calculations of the ratio of angle- of -attack changes 
resulting from torsioij to those resulting from bending; these ratios 
can be used to apply corrections to the results obtained through con- 
sideration of bending alone. The effects of the camber produced by 
both torsion and bending were neglected as these effects are relatively 
small. 

Inertia loadings are not considered in the analytical development 
and, therefore, the results apply only to airplanes which have a large 
percentage of their gross weight concentrated near the midspan. The 
alleviating effects of wing weight are considered briefly in the sec- 
tion entitled "Results and Discussion." 

In order to compute bending deflections under a given load ( or 
stress), the swept wings were transformed to unswept wings as shown 
in figure 1. This transformation is not exact for the case of wings 
having tapered plan forms; however, the error is small except for com- 
binations of low plan- form taper ratio and low aspect ratio. The root 
restraint for the transformed wing does not represent exactly the con- 
ditions for an actual swept wing because a streamwise twist due to 
bending which occurs near the root of a swept wing does not occur for 
the transformed wing. An analysis based on the assumed root restraint 
is believed adequate for the following reasons: 

(1) For many combinations of aspect ratio and sweep, the root twist 
is small compared to the angle -of -attack changes due to bending. 

(2) There is a wide variation in the root-restraint conditions in 
actual airplane designs. 

(3) Some experimental evidence exists which indicates that the 
bending deflections obtained for the transformed wing closely approximate 
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those for the swept wing when the effective root' is chosen in the manner 
used herein. 

Inasmuch as a direct algebraic solution for the aeroelastic effects 
considered herein was desired, it was necessary to express the spanwise 
load distribution responsible for the distortion in algebraic form. 

The feasibility of such a procedure was therefore investigated, and the 
results are presented in the appendix. These results indicate that the 
bending distortion of a wing supporting a given total load is fairly 
insensitive to the assumed form of the spanwise loading; in fact, the 
load shapes which result in approximately the same distortion are suf- 
ficiently broad to encompass the loadings to be expected on wings 
operating at widely differing flight conditions and having widely differing 
geometric parameters such as aspect ratio, sweepback angle, and taper 
ratio, as well as fairly large differences in degree of flexibility. 

On the basis of the results presented in the appendix, use of a 
representative average type of spanwise load distribution for purposes 
of computing the angle-of-attack changes resulting from bending deflec- 
tions appears to have at least limited application. The appendix also 
indicates that a uniformly distributed aerodynamic loading meets the 
requirements of a representative average type, and a uniform load modified 
by a tip correction was assumed for use in computing the angle-of-attack 
changes. Since the procedure used herein for examining aeroelastic effects 
primarily requires a knowledge of the difference between the loading on 
the rigid and flexible wings and since this difference was determined 
directly from the angle-of-attack changes due to the bending deflections, 
the exact details of the absolute loadings on either the rigid or flexible 
wing were not considered. The magnitude of the assumed uniform loading 
was determined by the operating lift coefficient of the distorted wing. 

This approach differs from the first step in the usual iterative or 
relaxation procedures (see ref. 2, for example), wherein the load on the 
undistorted wing is used to compute gin initial estimation of the distor- 
tion. The present approach affords a more accurate estimation of the 
distortion than is obtained from this first iterative step in that the 
magnitude of the loading assumed is correct and only its distribution is 
somewhat arbitrary. The angle-of-attack changes associated with the 
assumed loading were used to establish the difference in load between the 
rigid and flexible wings at any given section by assuming tnat this differ- 
ence was proportional to the product of the change in local angle of attack 
and the lift-curve slope of the rigid wing. 


METHOD OP AMIZfSIS 


Effects of wing bending .- The method of analysis involves relating 
the wing stress per g normal acceleration, the loss in wing lift due to 
bending, and/or the wing aerodynamic-center shift due to bending to the 
wing structural weight (these relations being a function of the wing 
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external geometry, structural material, and the flight condition) . For 
the transformed wing (see fig. l), the stress at any spanwise station is 
approximated hy the usual formula: 


a 


MrT’c* 

21 


(1) 


where 


■=' = c 'rf- U - X, i3 

T* = T' r [l - (1 - r)£j 
I = I 0 (T' ) 3 (c' )^(6r - 12r 1 2 + 8r3) 


p s' 

M B = J PC’ |(l - y‘) dl 

y 

C- S = c ’r[l - (1- 

As may he seen from the foregoing expression for thickness ratio, taper 
in thickness ratio was assumed to have a linear variation along the span. 

Performance of the indicated substitutions in equation (l) together 
with the following relations 

s' = A(1 + ,l) 
c'rT’r ^ T r cos 7 


T 

x r 


Eav_ 

1 + T 



results in an expression which relates the wing stress and loading to 
nondimens ional geometric parameters of the wing 


%/S 1 / A Y 
% I 0 \T a v cos y) 


128 


1 - (1 - \)£ 3 (6r - 12r 2 + 8r3) E - (l - t )?! 2 

S \ s 


(1 + t) 2 (1 + 1)2 


1 + 21 
~ 5 — 


1 + X. 


? - M - H 2 ® 


( 2 ) 
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For a given value of the left-hand member, hereinafter referred to as 
the "stress parameter," and for given taper ratios, equation (2) deter- 
mines the spanwise variation of the ratio of skin thickness to local 
airfoil-section thickness. Typical variations of a parameter propor- 
tional to skin thickness are presented in figure 2 for various plan- 
form and thickness taper ratios. The comparison applies to wings having 
the same area, sweep, aspect ratio, and average section thickness and 
having the same extreme-fiber stress under a given uniform loading. 

When the variation of skin thickness along the span is known, the wing 
structural weight for a given structural material can be determined from 
the relation 



(3) 


where 


a = 2a 0 T‘ (c* ) 2 r 


Performance of the indicated substitution in equation (3) together 
with the following relation 


(1 + \) d 

results in an expression for wing structural weight as a function of the 
variation of skin thickness along the span 



HA 1 / 2 

s3/ 2 wa Av 


16 

(1 + t)(1 + X) 2 





W 

The left-hand member of equation (4) is hereinafter referred to as the 
"weight parameter." Variations along the span of a parameter propor- 
tional to incremental wing weight are shown in figure 3 for various plan- 
form and thickness taper ratios. These wing-weight variations corre- 
spond to the skin-thickness variations of figure 2. 


Equations (2) and (^) can be combined graphically so that an evalua- 
tion of the weight parameter as a function of the stress parameter for 
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various values of the taper ratios X and t can "be made. This rela- 
tionship may he expressed as 

W A 1 / 2 
S 3 / 2 wa o-*-av 



,X,r 


( 5 ) 


The re la tionship between these two parameters is shown in figure 4 for 
various values of the taper ratios X and t . 

For the assumed wing structure and within the limitations of this 
analysis, the curves of figure 4 can he used to determine the wing 
structural weight for a given stress per g normal acceleration pro- 
vided the wing structural material, wing external geometry, and the 
airplane wing loading are known. Conversely, if the wing structural 
weight is specified, the corresponding stress per g normal accelera- 
tion may he determined. 

The stress parameter can easily he converted to parametric expres- 
sions for the aeroelastic effects considered herein. The ratio — 

a/n 

which appears in the stress parameter (see eq. 2) can he written in 
terms of the loss in wing lift or the wing aerodynamic-center shift due 
to bending by evaluating the angle -of -attack changes along the span 
resulting from the bending and then integrating the incremental loss 
in lift at each section and the associated incremental pitching moments 
that are thus produced about the aerodynamic center of the rigid wing. 

As stated in the foregoing development, the spanwise variations 
in skin thickness and wing structural weight were chosen so that a 
constant spanwise stress would occur in the extreme fiber under a uni- 
formly distributed load. With this type of skin- thickness distribution 
the stress along the span under an actual loading would be expected to 
be reasonably constant except in the vicinity of the tip where a gradual 
decrease in stress would occur due in part to the usual decrease in 
loading as the tip is approached. In addition, the stipulation of a 
constant stress out to the tip under any loading dictates that the skin 
thickness go to zero at the tip. At least a small finite skin thickness 
would exist at the tip for the practical case and, although such a modi- 
fication would not significantly affect the wing structural weight cal- 
culated by the preceding method, it would result in a decrease in stress 
to zero as the wing tip is approached. In order to account for these 
factors in the computations of the change in angle of attack due to wing 
bending, the wing stress was assumed constant to 0.8 semispan but was 
assumed to have - a linear variation from 0.8 semispan to zero at the tip. 
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When these assumptions are applied, the angle-of-attack change 
along the span may he obtained from the following relations: 


dz' 


and 


dzf 

dy 7 


Ax-d — sin 7 

dy' 


’ y ' r r/a ' 1 r 4 y 


Iq ET'c' 


E 


'0 


c'T* s'- 


For 0 < 2— < 0.8, 

a 1 ’ 


V = 


and 


( 6 ) 


dz' _ 2 as' 
dy' ET'j.c'j.oq 


y’A 


1 - (1 - t)4|[i - (1 - \) 
_ 8 L 


y' 


Integration of the preceding expression, substitution into equation ( 6 ), 
and subsequent rearrangement of the results yields 


^■B ET av 1 + r 1 + X 


oA tan 7 


HHSf-E 


(1 - X)J 


■ - !og e 1 - ( 1 - T 


(7a) 


For 0.8 < — < 1.0, 

o * 
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and 



Integration of the preceding expression, then substitution into equa- 
tion (6), and subsequent rearrangement gives 


oA tan 7 



0.8(1 - T )j - 



0 . 8(1 - 


V)] + 



(1 - 




lOg jl - (1 - T )| J- 


(7b) 


The parameter 


^B^av 
oA tan 7 


is plotted in figure 5 for various plan- 


form and thickness taper ratios. This parameter gives the variation 
along the span of the change in angle of attack produced by wing 
bending. A point worth noting is that the comparison of plan- form and 
thickness taper ratios shown in figure 5 is for wings having the same 
stress rather than the same structural weight. On this basis, the wings 
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become progressively lighter as the plan-form or the thickness taper 
ratio is decreased. 

The incremental lift associated with the angle -of -attack changes 
due to bending may be written in the following form: 


AL = 




The lift-curve slope 



associated with angle-of-attack changes 


due to bending is assumed to be the same as the lift-curve slope for 

the rigid wing in the region 0 < ^ < 0.8 but has a parabolic varia- 

s 

tion to zero at the tip in the region 0.8 < ~< 1.0 in order to 

s 

account for the usual decrease in lift as the tip is approached. 


For 0<Z< 0.8. 

s ’ 


( Ci <i)b = ( C Or 


For 0.8 < £ < 1.0, 
s ^ 


( c 0b- 5 W^ 
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When this assumption and the previously developed expressions for Ztog 
and c are used, the incremental load takes the following form: 



The incremental lift may he expressed in coefficient form and in terms 
of the nondimens ional geometric parameters of the wing through use of 
the following relations: 


AL _ AL = 

L " pS = C L 
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and 


sc 


r _ 


1 + X 


The substitution of these relationships into equation (8) and subsequent 
rearrangement gives 





5 



0.8(1 


(4 + _f_) log |T - 0.8(1 - T j] - 


x 0 + ( 5 + rrr) log E - (1 - x) l 




( 9 ) 


The expressions within the integrals on the right-hand side of 
equation (9) are proportional to the variation along the span of the 
change in lift due to wing bending. Typical spanwise variations of the 
change in lift produced by wing bending are plotted in figure 6 for 
various plan-form and thickness taper ratios. Note again that the com- 
parison of taper ratios is for wings having the same stress rather than 
the same structural weight. The sum of the integrals on the right-hand 
side of equation (9), which is used herein in determining the change in 
lift, is defined as f]_(X, t). The values of this function are plotted 
against plan-form taper ratio in figure 7 for r = 1.0 and t = 0.5. 
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Solution of equation (9) for the ratio of average loading to stress 
yields 


p _ w g/ s 

0 a/n 


tan 


7 


qk 



ET 


av 


Substitution of this expression into equation (5) relates a parameter, 
defined herein as the "incremental lift parameter," to the weight param- 
eter, that is. 


W 


. 1/2 


s 3/2 wa 0 T, 


= f 


av 


ac l /c l ei 


o^av) ( ^ 


tan 7 
® cos ^7 




(10) 


This expression is in the same form as the relation between the stress 
parameter and the weight parameter. Therefore, the incremental lift 
parameter may also be considered the abscissa of the plots shown in 
figure Within the limitations of the analysis, the curves of fig- 
ure 4 therefore can be used to determine the wing structural weight for 
a given loss in lift due to wing bending and vice versa. 

The lift-curve slope of the flexible wing may be expressed in terms 
of (see eq, 9) and the lift- curve slope of the rigid wing. For 

the flexible wing 


^ = ( C Ia )f 


For the flexible wing the lift equals the lift of the rigid wing plus 
the change in lift due to bending (negative for sweptback wings at posi- 
tive angle of attack) 


Lp — Lj^ + AT r — 




+ AL 


The preceding two equations can be combined and rearranged to give 



i C I<x)f _ Iff _ 1 

C I<x)r Iff - AL 1 _ AC^ 

C L 



( 11 ) 


Note that 


in itself is a function of 


(See eq. (10).) 
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The incremental pitching moment associated with the angle -of -attack 
changes due to tending may be written in the following form: 


AM 


= 2qs £ K)b ^B cd d ¥ 


where (see fig. l) 


d = y tan A - y tan A 


( 12 ) 


In spanwise coordinates 


y _ 1 1 + 2X 

T “ 3 ITT 


and 


. 1 1 + 
d = s 

\3 1 + 


— - — ) tan A 
X B 


When this expression for d and the expressions previously presented 
for Attg, (CLa^B* and c are substituted into the pitching-moment 

equation (eq. 12), the form of the equation is 


m 




ET a 


^(Tlhrf - |E - 11 - 


log [I - (1 - T 


>9} * * * ’ £° & hMe - <* - -a k - 5 © 2 • 3 {(“ 

log [I - 0.8(1 - t)] _ (u. + iog[i - 0.8(1 - xj] + ^5 - ttt) log £ - (1 - *■)£] 


( 5+ 5 1 

•) lOgjl - (1 - t)Z 

[ idjf) 

V t - 1 

/ L_ B - 

!i 7 


(13) 
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In order to express equation (13) in coefficient form and in terms of 
nondimens ional geometric, parameters of the wing, the following relations 
are used: 

/ 

m _ m = AC^qSc _ ACmg 

L pS C^qS Cl 

— — = 2. A i + ^ 

Sc ° 1 + X + X 2 

When these relationships are substituted into equation ( 13 ) and ‘the 
resulting equation is rearranged, the following expression is obtained: 


AC m P E ^av 1 _ r°* 8 _3_ (1 + X) 2 1 + T /l 1 + 2X 

Cp ct 5 A 2 (C I ^ R tan 7 tan A '-'0 1 ^1+X+X 2t-X '\ 31+X 



1 - (1 





15 (l + X) 2 1 + r /l 1 + 2X 

16 1 + X + X 2 T - X\3 1 + X 


1 - (1 




8 Z - 

s 





( 1 *) 


The expressions within the integrals on the right-hand side of 
equation (l4) are proportional to the variation along the span of the 
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incremental pitching moment due to wing bending. Typical spanwise 
variations of the pitching moments produced by wing bending are plotted 
in figure 8 for various plan-form and thickness taper ratios. The 
effect of taper ratio shown applies to wings of constant stress. The 
sum of the integrals on the right-hand side of equation (1*0., which is 
used herein to determine the pitching moment produced by wing bending, 
is defined as f 2 (l, T )- The values of this function are plotted in 
figure 7 against plan- form taper ratio. 

Solution of equation (1*0 for the ratio of incremental loading to 
stress yields 


a AC^/Cl E T av 2 ,T 


Substitution of this expression into equation (5) relates a parameter 
defined herein as the "incremental moment parameter" to the weight 
parameter; that is. 


W A 1 / 2 
S 3/2 wa 0 T av 


= f 


ACjh/Cl EI 0 T ay 3 


N 


tan 7 tan A f 2 (k, t) 


B 


p 

cos^y 


, X,T 


(15) 


This expression is in the same form as the relation between the stress 
and the weight parameter; therefore, the incremental moment parameter 
may also be considered as the abscissa of the plots shown in figure *J-. 
The factor £C m | Cp of the incremental moment parameter is the shift 

in the aerodynamic center of the wing (expressed in chord lengths) due 
to bending; that is. 



Within the limitations of this analysis the curves of figure k therefore 
can be used to determine the wing structural weight for a given allowable 
aerodynamic- center shift and vice versa.. 

Because the aerodynamic-center shift in terms of actual distance 
will vary with aspect ratio when the shift is expressed in chord lengths, 
for some purposes a better comparison is afforded when the shift is 



22 


NACA TN 2901 


expressed in terms of (based on over-all size of the airplane). 

The relationship between Ax/c and Ax J^S may be developed from the 
expression 

- _ ij- 1 + x + x 2 IF 

3 (1 + X) 2 » A 


so that 


Ax _ Ac _1_ 4 1 + X + X 2 

ifs c tfA 3 (1 + X) 2 


The function 


4 1 + x + x 2 

3 (1 + X) 2 


is defined herein as 


against plan-form taper ratio in figure J. 


f 3 U) 


and is plotted 


The location of the center of pressure of the change in load pro- 
duced by wing bending is another- parameter useful for the study of the 
effects of wing bending on the longitudinal stability of a complete 
airplane. The streamwise location of this center of pressure expressed 
in chord lengths from the aerodynamic center of the rigid wing is 
ACjjJaC^ which from a combination of equations (9) and (l4) can be 
expressed as follows: 


ACL, 

S = A tan A 


f 2 (X,r) 


The ratio 


f2(\>T) 

fl(Vr) 


is plotted in figure 7 as a function of plan- form 


taper ratio for thickness taper ratios of 1.0 and 0.5. When the aero- 
dynamic center of the rigid wing is assumed to be located at one-quarter 
of its mean aerodynamic chord, the expression for the center of pressure 
of the change in load is 



-0.25 + A tan A 


f 2 (X,r) 

fl(Vr) 


c 


(16) 
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Effects of wing torsion .- The variation of the torsional deflection 
of a wing along its span may he expressed in the following form: 


ae = % = % = % a 2 

dy' GJ G — I ® T' c* — 

I I 


(IT) 


The twisting moment may he expressed as 


% = / p(c’ ) e dy' 

Jy» 


= (c* )‘ 


p s ' | — T1 2 

L (1_ ay ’ 


(18) 


When a constant loading and a constant value of e along the span are 
assumed, the integration of equation (l8) yields 


Mj, = pe(c , r ) 2 s' 1 


1 + 1 + X* 


y 

s 


(1 . x) (zf . ^(i) 3 


(19) 


For a constant loading condition the bending moment may he expressed 
as (see expression for Mg following eq. l): 


Mg = pc , r (s t ) 2 


1+21 _ l+l y + 1 (jf 1 - l /y\ 3 

6 2 s 2 \B/ 6 \s) 


( 20 ) 


Dividing equation (19) hy equation (20) to obtain the ratio of the 
twisting moment to bending moment gives the following expression 
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(21) 
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In the analysis of the effects of torsion only wings with constant thick- 
ness ratio (t = l) are considered. When this assumption of constant 
thickness ratio is made, equation (21) is substituted into equation (17) 
and the result is rearranged and integrated, the torsional deflection 
may he expressed as 


0TG j cos 7 
2 ae 


py/ s 


uo 


1 + X + X 2 


(1 - 

x) 2 (y\ 3 


L 3 

3 




’ + 2X 1 + \ y ( 
6 2 s 
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eg 


y 
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For the case of zero taper ratio, equation (22) simplifies to 


6TG- j cos 7 
2 ae 




y 

s 


- -2 l°Be(l - i) 


(23) 


The torsional deflection is related to the angle of attack due to twist 


AOm 

e = — — 

COS 7 


Substituting equation ( 2k ) into equation (23) gives 




(24) 


(25) 


The angle-of-attack change due to bending, Axg for 0< —< 0.8, 

is given by equation (7(a)). Substituting X = 0 into equation (7(a)) 
gives the following equation: 


— \ iog e ii 


* 



aA tan 7 


(26) 
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Dividing equation (25) "by equation (26) and solving for the ratio of 
angle-of -attach change due to torsion to the angle-of-attack change due 
to bending yields 


zkirp 8Ee 

AG ^ tan 7 

Similar computations for X = 1.0 and r = 1.0 give 

= k lQ ge(l - g)Ee 

AG — — tan 7 
I s 


(27) 


(28) 


RESULTS AND DISCUSSION 


Comparison of effects of wing bending and wing torsion .- An exami- 
nation of equations (27) and (28) shows that the ratio of the angle-of- 
attack changes due to torsion to those due to bending are inversely 
proportional to the aspect ratio and the tangent of the angle of sweep- 
back. The over-all variation with sweepback would be approximately in 
inverse proportion to the angle of sweepback for small angles, but the 
angle-of-attack ratio would decrease rapidly at large angles of sweep- 
back. These equations also indicate that the angle-of-attack change 
due to torsional deflections always opposes those due to bending when- 
ever the center of pressure is ahead of the elastic axis. The angle- 
of-attack ratio for a plan-form taper ratio of zero Is constant along 
the span; whereas for a taper ratio of unity the ratio has a variation 
along the span from a value one-half that for zero taper ratio at mid- 
span to infinity at the tip. The presence of this value of infinity 
at the tip results from the fact that the stress was assumed to remain 
constant out to the tip when the torsional deflections were considered. 
(See eq. 23.) As was pointed out previously in the practical case the 
bending stress would decrease near the tip to zero at the tip. Because 
of this consideration, the use of equation (28) should be limited to 
regions inboard of the 0.8 semispan as was done when bending deflection 
alone was considered. At the 0.8-semispan station the value of the 
ratio of change in angle of attack would be approximately the same for 
the unit taper ratio as for zero taper ratio. 

In order to examine the actual ratios of angle-of-attack changes 
due to bending and torsion, it is necessary to know the values of the 
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ratios g/e a-n^ j/l. The ratio g/e is about the same for both steel 
and duralumin and has a value of about 0.37* The ratio j/l may have 
rather large variations and still maintain adequate strength in the wing. 
In fact, it is theoretically possible to adjust the torsional moment of 
inertia of a wing to obtain complete compensation of aeroelastic effects 
at a given flight condition (isoclinicism) . In most cases, however, the 
torsional moment of inertia is dictated by other considerations such as 
those of flutter or control effectiveness. 

In order to investigate the ratio of j/l inherent in an airfoil 
section having the type of structure assumed in the present report, the 
torsional moment of inertia of a symmetric parabolic -arc section was 
calculated from the relation 



(29) 


where 



is a line integral taken around the mean periphery of the 


torsion box consisting of the upper and lower skins and the shear webs. 


The contribution of the areas forward of the front shear web and 
back of the rear shear web to a was neglected because it was relatively 
small and the contribution of the shear webs to the line integral (thin 
airfoil sections) was also neglected. For these conditions equation (29) 
may be written for a symmetric parabolic-arc section as follows: 


where 



2(a) 2 r 



(30) 


f(x’) = iffi'c* — (l - — ) 
c* V c v 
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a = c 


(1 - r) 


>1- 4 

c’ 


ik/c 


f(x’ ) d 4 
s ’ 


= | T* (1 - r 


„ c . )2 |r. 6 (i) 2 +4 (i.g 




When the indicated substitutions into equation (30) are made and k is 
assumed to be equal to 0.2 (shear webs located at 20- and 80-percent 
chord ) , the value of the torsional moment of inertia becomes 


J = 0.803(T t ) 3 (c ' )\r - 2r 2 + r3) (31) 


Using the expression for I previously presented (see expression fol- 
lowing eq. l) and substituting for Iq its value for a symmetric 
parabolic-arc section (O.O386) results in the following expression for 
the ratio j/l 


- = 3.52 — =-^- 

1 1 - 2r + | r 2 

This ratio has some variation with skin-thickness ratio but the varia- 
tion is slight except when the section approaches solidity; the value 
oi j/l for the assumed section is about 3.5» This value is reasonably 
typical of values existing for actual wing sections. 

When the value of j/l (3-5) and the value of G/E previously 
discussed (0.37) are used, equations (27) and (28) may be written 


(32) 
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for X = 0, 


1 


Amp 6.2e 

Axjj A tan 7 


(33) 


and for X = 1.0, 

3.1 l°8e(l - f) e 

Ax B y/s A tan 7 


The- values of e for a sweptback wing are not constant along the span, 
but the average value probably ranges from a maximum positive value of 
about 0.25 at subsonic speeds to small positive or negative values at 
supersonic speeds. By assuming a conservative value for e of 0.25, 
the ratio of angle of attack due to torsion to that due to bending has 
been plotted in figure 9 as a function of taper ratio and sweep angle. 

This figure shows that, for an angle of sweepback of the elastic axis 
of 45° and an aspect ratio of 4, the ratio of angle-of-attack change due 
to torsion to that due to bending would be about -0. 38 for the zero- 
taper-ratio case. Reducing the aspect ratio to 2 results in this value 
being doubled and reducing 1 * the sweepback angle to 22.5° has a somewhat 
greater effect. Conversely, doubling the aspect ratio to 8 would halve 
this value, whereas increasing the sweepback angle would reduce this 
value much more rapidly. As mentioned previously, the values of the 
ratio for the constant-chord case are one-half the value for the zero- 
taper-ratio case at the root and about the same at the 0.8-semispan sta- 
tion. The variations presented in figure 9 indicate that, except for 
small sweepback angles (the order of 30°) or small aspect ratios (the 
order of 2) or particularly a combination of .the two, the angle-of-attack 
changes produced by torsional deflections are appreciably less than those 
produced by bending deflections. In view of this result the remainder 
of this paper is concerned primarily with the effects of wing bending. 

If the ratio of torsional moment of inertia to bending moment of 
inertia can be reduced from the value of 3*5 assumed in this analysis, 
the compensating effect of torsion would be increased. The distance 
between the aerodynamic center and the elastic axis also affects the 
amount of torsional compensation, and, in many cases, this distance will 
be less than the value of 0.25 chord used in the present analysis. A 
reduction in this value tends to reduce the compensating effect of torsion. 

General effects of wing bending .- The effect of wing structural 

parameters , wing geometric parameters, and flight condition on the rela- 
tion between the structural weight of a wing and aeroelastic effects 
dependent on wing bending stiffness may be examined through use of 
figure 4. 
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The structural-weight parameter, the ordinate of figure 4, indi- 
cates that the weight of a wing varies as the cube of a linear dimen- 
sion (or as S 3/2). This variation is the type obtained with geometri- 
cally similar structures. The structural-weight variation with size 
for actual wings is found to be slightly less than the cubic relation 
because the efficiency of wing structures increase with increase in 
size, at least up to a certain point. These actual variations, however, 
are not reduced from the cubic variation a sufficient amount to prevent 
the weight of wings designed to a given stress from becoming prohibitive 
as their size is increased. The reason that the weight becomes prohibi- 
tive is that the lift-producing capabilities of a wing increase only as 
the square of a linear dimension. In order to circumvent this diffi- 
culty when the structural design of wings is dictated primarily by 
stress considerations, designers tend to reduce the allowable load fac- 
tor as airplanes become larger; thus, the tendency for the weight to 
increase is retarded. When factors such as aerodynamic-center shift 
due to wing bending are a consideration, such an alternative is not 
afforded the designer because there is no reason to expect that the 
aerodynamic -center shift may be reduced significantly as the airplane . 
size increases. 

The stress parameter shown in figure 4 indicates that the allowable 
stress relative to the density of the structural material is a primary 
consideration in selection of the material when the wing structural 
design is dictated by strength alone. When stiffness governs a design, 
however, the incremental moment and lift parameters shown in figure 4 
indicate that the modulus of elasticity of the structural material rela- 
tive to its density also would be a primary consideration. 

The incremental lift and moment parameters show that, as far as the 
effect of flight condition on aeroelastic effects is concerned, dynamic 
pressure is the primary variable. Mach number has only an indirect 
effect on aeroelastic phenomena. Changes in Mach number produce changes 
in the Incremental lift or moment parameters through the effect of Mach 
number on the lift-curve slope of the wing. The relation between the 
lift-curve slope and the Mach number, aspect ratio, and sweep angle 
assumed for use in the examples presented herein is given by the 
expression 

c L a = ~ p ■■■ 2 * 

\l - M^cos^A + 2 
cos A A 

where M is the Mach number. This expression is that given in refer- 
ence 5 with the addition of the Prandtl-Glauert correction for compres- 
sibility which has been applied to the two-dimensional lift-curve slope. 
The formula gives values of lift-curve slope in reasonable agreement with 
those determined experimentally . A more rigorous determination of lift- 
curve slope from a simplified lifting-surface theory can be obtained 
through use of reference 6. 
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The stress, incremental lift, and incremental moment parameters 
reveal that the strength and stiffness characteristics of a wing are 
strongly affected by its external geometry. Increase in aspect ratio, 
sweep angle, or decrease in section thickness ratio caused rapid (high 
order) increases in stress and even more rapid increases in aeroelastic 
effects. In addition to the .direct effects of aspect ratio on the param- 
eters presented in figure the lift-curve slope is affected by aspect 
ratio. The effect of an increase in .aspect ratio on the wing lift- 
curve slope is in a direction to increase aeroelastic effects, and 
the effect of an increase in sweep angle on the wing lift-curve slope 
is in a direction to decrease aeroelastic effects. Although the 
area and moment of inertia of the airfoil section appears in the param- 
eter presented in figure b, the effects of the geometry of the airfoil 
section are generally small because of the limited variation of airfoil 
shapes which provide satisfactory aerodynamic characteristics. The 
variations in nondimens ional moment of inertia I Q of airfoil sections 

are larger than -variations in nondimens ional area ag but are much less 

significant than the effects of variations in the thickness ratio of 
these sections or the plan-form geometry of the wing. 


Effects of wing external geometry .- In order to provide a better 
illustration of the manner in which various geometric properties of a 
wing affect its strength and its aeroelastic characteristics, the charts 
presented in figures ^ and 7 have been used to calculate the Btress 
per g normal acceleration, the ratio between lift-curve slopes of the 
flexible and rigid wing, and the aerodynamic-center shift due to bending 
(in terms of both chord lengths and for a series of duralumin wings 

flying at a Mach number of 0.9 at an altitude of 30,000 feet. All the 
wings have symmetric parabolic-arc sections, a structural weight typical 

of current practice (w/S^/^ = 0 . 2 ), a plan-form taper ratio of 0.5, and 
no variation in section thickness ratio along the span. Wings for which 
the change in. lift due to bending was greater than half the lift remaining 


on the distorted wing 


( c O 


R 


< ^ were not included in these calculations 


because limitations imposed by the assumptions Used in the present analysis 
do not warrant such calculations. The preceding value was selected as a 
limit because a study indicated that, for the cases where losses in lift 
are larger than this value, the variation in load distribution from the 
assumed uniform type would, in some cases, be more extreme than the ones 
investigated in the appendix. As previously stated in "Assumptions and 
Limitations of Analyses" these investigated variations from the assumed 
load distribution have only a secondary effect on the predicted change 
in lift between the flexible and the rigid wing. When not limi ted by 
the foregoing restriction, calculations were made for aspect ratios up 
to 10 and sweep angles up to 65°. 
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The results of the calculations for this family of wings are sum- 
marized in figure 10 where the computed stress and aeroelastic param- 
eters are plotted against sweepback angle for' various aspect ratios and 
section thickness ratios- The stress per g is plotted in the form 
shown because the stress is determined by the magnitude of the loading 
Wgn/S rather than by the acceleration alone. The effect of such varia- 
tions is shown subsequently. 

For wing structural weights typical of current practice, the varia- 
tions shown in figure 10 indicate that for a very low thickness ratio 
(T = 0.04) only low aspect ratios (the order of 2) are usable for swept 
wings, and then both stress and aeroelastic effects will be important 
at large sweep angles (the order of 60°). Moderate thickness ratios 
(T = 0.08) apparently can be used in conjunction with moderate aspect 
ratios (up to roughly 6) at small angles of sweepback (about 30° )> but 
the usable range of aspect ratios decreases rapidly as the angle of sweep 
increases. With a thickness ratio more typical of current practice 
(T = 0.12) fairly high aspect ratios (of the order of 8) are possible 
with small amounts of sweepback (about 30°) but low aspect ratios appear 
necessary if large amounts of sweep are desired. 

As mentioned previously, the charts presented in figure 4 were 
obtained from considerations of wing bending alone; therefore, the 
variations presented in figure 10 and subsequent figures do not include 
the effects of torsion. The extent to which torsional distortion will 
modify the trends shown may be judged by reference to figure 9 or to 
equations (2J) and (28), which define the ratio of the change in angle 
of attack due to torsion to the change in angle of attack due to bending. 

The equations rather than figure 9 should be used when the values of 
j/l or e differ significantly from those on which figure 9 Is based. 
Correlation of the trends presented in figure 10 and the effects indi- 
cated by figure 9 leads to the conclusion that, whenever the loss in lift 
or shift in wing aerodynamic center produced by wing bending becomes 
large, the percentage alleviation afforded by torsional distortion is 
usually not great. For example, if wings with sweep angles greater 
than 30° are considered, wings with combinations of aspect ratio, thick- 
ness ratio, and sweep which are shown by figure 10 to result in aerodynamic- 
center shifts of 10 percent or greater due to wing bending would .in no 
case have more than half of the shift compensated by torsional distortion 
and, in most cases, the compensation would be much less than this value. 

It should be remembered that the stress and aeroelastic effects could 
also be alleviated somewhat frpm those presented through use of a higher 
taper ratio or wing structural weight. 

As indicated by the variations shown in figure 10, increases in 
either aspect ratio or sweep angle produce much larger percentage 
increases in aerodynamic -center shift due to bending than in stress. 
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This result shows that flexibility (in terms of its effect on longi- 
tudinal stability) tends to become more significant relative to stress 
at high aspect ratios and sweep angles. 

Results are presented in figure 11 for a family of wings similar 
to that presented in figure 10 but the thickness ratio of which was 
held constant at 0.08 and the plan-form taper ratio of which was varied. 
The computed stress and aeroelastic parameters are plotted in figure 11 
against aspect ratio for various values of sweep angle and plan-form 
taper ratio. These results indicate that for constant-chord wings 
aeroelastic effects are appreciable for all but fairly low aspect 
ratios. Decreasing the taper ratio results in a significant increase 
in the combinations of aspect ratio and sweepback which will not exceed 
given values of the aeroelastic parameters. 

The mann er in which plan- form taper ratio affects the relative 
importance of stress and the two aeroelastic effects considered herein 
may be investigated through reference to figure 7- The function f^(A,/r) 
exhibits a continuous increase with increase in taper ratio. Since this 
function appears in the incremental lift parameter but not in the stress 
parameters (see fig. 4) a progressively greater loss in lift due to wing 
bending relative to the stress per g normal acceleration is indicated. 
Correspondingly, the function f2(l/ r ) which appears only in the incre- 
mental moment parameter has a rapid increase up to a taper ratio of 0.4 
but little change thereafter. (See fig. 7*) This variation means that 
the aerodynamic -center shift due to wing bending becomes increasingly 
larger relative to the stress per g with increase in taper ratio at 
small taper ratios but there is little change with further increase in 
taper ratios. The variation of the ratio f2( A.,t )/f-^( A,,t ) indicates 

the relative importance of aerodynamic -center shift and loss in lift as 
aeroelastic effects. Thus, the magnitude of the aerodynamic-center 
shift relative to loss in lift increases with increase in taper ratio 
up to a taper ratio of about 0.25 but then decreases with further increase 
in taper ratio. The relative magnitude of these two aeroelastic phenomena 
may be important in establishing the over-all effect of aeroelastic ity on 
the longitudinal stability of an airplane with a horizontal tail. This 
point is discussed subsequently. 

The effect of a linear variation in thickness ratio along the span 
is indicated in figure 12 where the computed stress and aeroelastic 
parameters are plotted against angle of sweepback for values of r 
equal to 0.5 and. 1.0. The wings used in this example have the same 
section and structural weight, operated under the same flight condition, 
and used the same structural material as the previous examples. These 
wings have an aspect ratio of 4, a symmetric parabolic-arc section, and 
the average of the root and tip section thickness ratios was 0.08. 
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A linear reduction in thickness ratio to a value at the tip of 
one-half that at the root, as shown in figure 12, reduces the stress 
and aeroelastic effects by a small amount over those which occur for 
the wing with a constant thickness ratio. Whether advantage could he 
taken of such a modification would depend upon its effect on other 
contributing factors in the wing selection such as drag. Reference 
to the variations in f^(l,r), f2(l,T), and \,t) ff-^i X,t) in fig- 

ure 7 indicates that tapering the thickness ratio produced a somewhat 
smaller ratio of loss in lift to the stress per g normal acceleration 
than occurred for the constant-thickness -ratio case and a slightly 
larger ratio of aerodynamic-center shift due to bending to stress per g 
normal acceleration than occurred for the constant-thickness-ratio case. 
The shift in aerodynamic center relative to the loss in lift is larger 
than for the constant-thickness -ratio case. 

Effect of wing structural weight .- The effect of changing the wing 
structural weight is shown in figure 13. The duralumin wings assumed 
in this example had an aspect ratio of 4, a constant thickness ratio 
of 0.08, a plan-form taper ratio of 0.5, a symmetric parabolic-arc sec- 
tion, and operated at the same flight conditions as in previous examples. 
The magnitudes of the stress or aeroelastic effects are approximately in 
inverse proportion to the structural weight for the range of structural 
weights presented. For higher structural weights, these variations 
would be less because the wings approach a solid condition, and the 
increased structure afforded by the increase in weight must be added 
nearer the elastic axis where it does not contribute appreciably to the 
strength and stiffness of the wing. 

Effect of structural material .- The effect of the choice of steel 
or duralumin as a structural material is show in figure l4. The example 
wing and the flight condition are the same as’ those considered in the 

preceding section and the value of w/S^^ selected was 0.2. As may be 
seen from figure l4, very little choice can be made between these 
materials as far as aeroelastic considerations are concerned, although 
duralumin is usually considered superior on the stress -weight basis. 

This result with respect to aeroelastic effects is due to the fact that 
the increased density of steel over duralumin is offset by an almost 
proportional increase in modulus of elasticity. Actually, a superiority 
of steel wings over duralumin wings is indicated on the basis of aero- 
elastic considerations when the duralumin wing approaches solidity 
because the steel wing of the same weight would still be relatively 
hollow and therefore more efficient structurally. 

Effect of flight condition .- The effect of changes in-flight condi- 
tion is shown in figure 15. The duralumin wings assumed in this example 
had an aspect ratio of h, a constant thickness ratio of (3.08, a plan- 
form taper ratio of 0.5, a symmetric parabolic-arc section, and a 
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structural weight corresponding to a value of W/S 3/2 of 0.2. The com- 
puted stress and aeroelastic parameters are plotted against angle of 
sweepback. The curves shown in figure 15 are for three flight condi- 
tions. The basic curve is for the condition used in previous examples, 
a Mach number of 0.90 at an altitude of 30,000 feet. A second curve 
shows the effect of a reduction in altitude to sea level while main- 
taining the same airspeed, and a third curve shows the effect of a 
speed change at an altitude of 30,000 feet corresponding to a reduction 
in Mach number from 0.90 to 0.75. The upper plot on figure 15 indicates 
that there are no effects on the stress per g normal acceleration 
because in the method of analysis no account was taken of the effect on 
the stress of the inboard shift of the center of pressure due to aero- 
elasticity. In most cases, this effect would be small. The large 
increase in aeroelastic effects due to a decrease in altitude results 
from the change in dynamic pressure associated with the increased 
density. The decrease in Mach number had a slight effect for the sea- 
level condition. This effect is manifested by a reduction in lift-curve 
slope and is relieving in nature, as indicated by the equation for Cj^ 

presented previously. The figure also indicates a relatively large 
decrease in aeroelastic effects for a small decrease in airspeed. Again 
this reduction primarily reflects a decrease in dynamic pressure. A 
secondary effect results from reduction in lift-curve slope through 
reduction in Mach number. 

Effect of inertia loading .- No detailed discussion can be presented 
as to the effects of the inertia of the wing on aeroelastic phenomena. 
These effects cannot be generalized in that they depend on the distri- 
bution of wing weight which can have extremely wide variations due, for 
example, to the presence of concentrated masses that do not contribute 
to the structural strength. The inertia of the wing, however, does 
produce a relieving effect in maneuvers which may be illustrated by 
assuming that the total wing weight is uniformly distributed along the 
span (the same assumption as for the aerodynamic loading). Under these 
conditions the inertia loading p^^ is given by 


Pi = - 



and the aerodynamic loading p a is given by 


Pa 


W^n 

S 
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The total loading then is 



For most airplanes having a fuselage and relatively thin wings W w / W q. 

is small compared to unity. This condition gave rise to the assumption 
used herein that only the aerodynamic loading (first term in the fore- 
going equation) he considered in the present analysis. The foregoing 
equation shows/ as is well-known, that the wing weight tends to reduce 
the loading which produces the aeroelastic distortion. In cases where 
the weight contained within the wings is large compared with the gross 
weight (the flying wing, for example), this relieving effect would he 
extremely beneficial. Advantage may he taken of corresponding effects 
of concentrated loads through judicious location of engine nacelles and 
external stores. 

Effects involving the complete airplane .- The previous discussion 
has been concerned with the effects of aeroelasticity on the wing itself. 
When the over-all effects to he expected on an airplane with a horizontal 
tail are considered, further discussion is appropriate. The results 
obtained herein for a wing can he applied to a horizontal tail as well, 
hut other effects on the complete airplane, such as bending of the fuse- 
lage or the induction effects caused by the downwash changes at the tail 
produced by wing distortion, cannot he analyzed by the present method. 
These latter effects are usually smaller than the effects caused by 
distortion of the wing and tail, however, and the discussion is thus 
confined to these two components. Effects of both horizontal-tail and 
fuselage flexibility on longitudinal control are considered in refer- 
ence 7 and similar effects on longitudinal stability are also considered 
briefly. 

The wing-aerodynamic -center shift caused by aeroelastic effects 
directly reflects a change in the longitudinal stability of a flying- 
wing airplane. The change in the longitudinal stability of an airplane 
with a horizontal tail is in turn related to a shift in the point for 
neutral stability with respect to angle of attack. This neutral point 
is not only affected by the wing-aerodynamic-center shift but also .by 
other factors. This condition may be seen from the following equation 
which expresses this shift 
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The pertinent assumptions made in deriving the foregoing expression are 
that the aerodynamic -center shift due to wing tending Aoc/c - is the same 
for the wing-body combination as for the wing alone, that the distance 
between the point for neutral stability and the tail center of pres- 
sure Z is the same for the rigid airplane and the flexible airplane, 
and that the variation of downwash angle at the tail with angle of attack 
de/da is not affected by the wing distortion. This equation shows that 
the neutral-point shift is dependent on a term involving the effect of 
flexibility on the lift-curve slopes of the wing and horizontal tail as 
well as a term involving the shift in aerodynamic center of the wing. 

For the case where ( C Ic 4 : ) F ^( C Ia i .) R 18 e< 3 ual to ( Cl <* w ) F 

(a condition which might conceivably be called equal flexibility of wing 
and tail), the foregoing equation shows that the shift in neutral point 
is the same as the shift in the aerodynamic center of the wing. For the 

case where is greater than 1^ < the tail 

effectively stiff er than the wing), the second term in the foregoing 
equation becomes negative and tends to compensate for the positive shift 
in wing aerodynamic center produced by aeroelastic distortion. Presum- 
ably, if such a condition exists, it should be possible to adjust such 
factors as S-^/S or Z/cT until the wing-aerodynamic-center shift is 

completely compensated. The physical interpretation of the condition is 
that the wing would be so positioned on the fuselage that the center of 
pressure of the loss in lift for the complete airplane caused by aero- 
elastic action would be at the point of neutral static stability for the 
rigid airplane, and therefore no shift in this point would occur because 
of flexibility. Generally, it would be more feasible to keep the tail 
length Z constant and to adjust the tail size to accomplish this end. 

If the horizontal tail is assumed to be rigid the tail size required 
for complete compensation of the aeroelastic effects considered can be 
obtained by equating the expression for the location of the center of 
pressure of the loss in wing lift previously derived (eq. (l6)) to the 
expression for the center-of-gravity position for neutral static sta- 
bility of the rigid airplane. With a flexible tail, a different 
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relation would apply, and the effect of tail flexibility is discussed 
in more detail subsequently. Equating the expressions as stated pre- 
viously gives 






de\qt ^t l 

day q S c 
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In order to express the various lengths appearing in this equation in. 
terms of a general length dime nsion rather than the chord, the following 
expression is substituted: 



Solving for the qrea ratio gives 
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When more or less typical values for the parameters in this equation 
are assumed, area ratios have been calculated for various aspect ratios 
and singles of sweepback. These area ratios are presented in figure 16. 
The assumed values of the parameters are as follows: 
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The value of 0.25 — + — is typical for a wing-fuselage combination 

rather than for a wing alone. The center of pressure of the loss in 
wing lift is also presented in 'figure 1 6 as a function of aspect ratio 
and sweephack angle. 

The results presented in figure 1 6 show that the tail size required 
for complete compensation of the wing aeroelastic effects on longitudinal 
stability increases somewhat with increase in aspect ratio. The tail 
size increases much more rapidly with increase in sweephack angle, par- 
ticularly at large sweephack angles. In general, the tail sizes required 
for this complete compensation would appear prohibitive except at small 
angles of sweephack. Usual tail-wing area ratios are of the order 
of 0.2 and the aerodynamic efficiency of an airplane is known to decrease 
as the size of the horizontal tail increases. Reference again to fig- 
ure 16 indicates that for most conditions the center of pressure of the 
loss in wing lift is further hack than the neutral -point locations for 
most airplanes. Location of the point for neutral stability with respect 
to angle of attack usually ranges from 30c to 50 c hack of the 'leading 
edge of the mean aerodynamic chord. It does not appear, however, that 
rearrangement of the airplane configuration to provide satisfactory static 
margins with extremely rearward center-of-gravity positions would he par- 
ticularly difficult from the standpoint of rearranging the weight dis- 
tribution; the large tail-size requirement would he the most serious 
limitation. The possibility for obtaining complete compensation appears 
even more remote when a flexible tail is considered. For any given wing 
sweephack angle and aspect ratio the tail size required for complete com- 
pensation increases rapidly with increased tail flexibility until, when 
the wing and tail are equally flexible, no compensation for the wing- 
aerodynamic -center shift is possible regardless of tail size. As would 
be expected if the horizontal tail is more flexible than the wing, the 
presence of the tail is detrimental from the aeroelastic standpoint and 
increase in tail size would result in increased effects of aeroelasticity 
on the longitudinal stability of the airplane. On the other hand, when- 
ever the horizontal tail is effectively stiffer than the wing, some degree 
of compensation is always obtained. 


CONCLUSIONS 


An analysis has been made of factors affecting the loss in lift and 
the shift in aerodynamic center of a wing produced by the wing bending 
under aerodynamic load. The analysis was applied to shell-wing structures 
in which the spanwise variation in skin thickness is of a character to 
give a constant spanwise stress under a uniformly distributed load. 
Conclusions obtained from the results of this analysis are as follows: 
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1. The geometric parameters the variations of which produce the 
largest changes in the aeroelastic effects considered herein are aspect 
ratio, sweepback angle, and section thickness ratio. Based on the 
weight of current wing structures, consideration of the effects of wing 
bending alone leads to the conclusion that for sweptback wings with low 
values of section thickness ratio (0.04) the aeroelastic effects are 
extreme for all but low aspect ratios (2) and even then the effects are 
large for large angles of sweepback (60°). The effects of wing flexibility 
on the longitudinal stability of an airplane become more important rela- 
tive to considerations of wing stress as the aspect ratio op sweepback 
angle of the wing is increased or the thickness ratio decreased. 

2. The effects of torsional flexibility on loss in wing lift and 
shift in aerodynamic center will usually tend to alleviate the effect of 
bending, but the alleviation afforded will not be large in general because 
the angle-of-attack changes due to torsion will usually be much smaller 
than those due to bending except for wings with low values of sweep, 
aspect ratio, or, in particular, a combination of the two. 

3. Reducing the plan-form, taper ratio reduces the aeroelastic effects 
considerably. The reduction in the shift in aerodynamic center resulting 
from decrease in taper ratio is about the same as the reduction in stress 
for all but very low taper ratios where the reduction in the shift of the 
aerodynamic center is much greater than the reduction in stress. The 
reduction in loss in lift through reduction in taper ratio is larger than 
the reduction in stress for all ranges of plan-form taper ratio. 

4. For wings having a constant extreme -fiber stress along the span, 
the alterations in aeroelastic effects produced by variations in section 
thickness ratio along the span and by practical variations in airfoil 
shape are small. 

5. Of the factors dependent on flight condition, dynamic pressure is 
the most important from the standpoint of aeroelastic effects. For the 
cases considered, Mach number effects exist chiefly in the manner in 
which the lift-producing effectiveness of the wing is altered. 

6. An increase in wing weight per unit area with increase in size 
appears to be more difficult to prevent when aeroelastic effects are an 
important consideration than when stress alone is a consideration. For 
the type of structure considered herein and for any given wing material 
and external geometry, the loss in lift and aerodynamic- center shift are 
approximately in inverse proportion to the structural weight except when 
the wing is nearly solid. Under the latter condition, the effects of 
changes in wing structural weight are reduced. 

7. On the basis of the aeroelastic effects considered, there 
appears to be little choice between steel and duralumin as a structural 
material except for conditions when the duralumin wing approaches solidity; 
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the steel wing of the same weight would still he relatively hollow and 
therefore superior because of its higher structural efficiency. 

8. In man euvers the aeroelastic effects produced by the aerodynamic 
loading are alleviated somewhat by the inertia of the wing. 

9. The effects of loss in wing lift and shift in its aerodynamic 
center on the static longitudinal stability of an airplane with a hori- 
zontal tail are dependent on the flexibility of the horizontal tail. If 
it is assumed that the fuselage flexibility is small and that the effect 
of wing flexibility on the variation with angle of attack of the downwash 
at the tail can be neglected, the following conclusions apply: 

(a) The shift in wing aerodynamic center is approximately equal 
to the shift in the cent er-of -gravity position for neutral stability 
with respect to angle of attack when the loss in lift of the wing and 
tail are equal. 

(b) The shift in wing aerodynamic center is usually somewhat greater 
than the shift in center-of-gravity position for neutral stability with 
respect to angle of attack when the loss in lift of the tail is less 
than that of the wing. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, August 15, 1952. 
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APPENDIX 

THE EFFECT OF SPANWISE DISTRIBUTION OF LOADING ON THE 
SPANWISE VARIATION OF ANGLE-OF -ATTACK 
CHANGE DUE TO BENDING 


In order to illustrate the effect of various distributions of span- 
wise loading on the angle of attack due to bending, two extreme loading 
conditions were investigated and the angle -of -attack variations were 
obtained and compared with those obtained for a constant loading condi- 
tion assumed in the present report. The conditions investigated are 
(l) the load varies linearly from root to tip, the tip-chord loading 
being 50 percent less than the root-chord loading, and (2) the load 
varies linearly, the tip-chord loading being 50 percent greater than the 
root -chord loading. Parameters applying to the first extreme-loading 
condition are denoted by the subscript 1; parameters applying to the 
second extreme-loading condition are denoted by the subscript 2; and 
the parameters applying to the constant loading condition have no sub- 
script. For these calculations a constant thickness ratio was assumed. 
In order to analyze the problem, the following relations were used to 
determine the ratios of the angle-of-attack changes along the span 
(see expression following eq. 6): 


Zkt B = - 


dz* 

dy' 


sin 7 


where 




The ratios of angle-of-attack changes are then found as 



(Al) 
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The stress ratios are proportional to the ratios of the bending moment 
and are 




% 

cr v t 
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Since the wings were designed for a constant spanwise stress for the 
constant -loading condition, equation (Al) may be expressed 



The moment at any spanwise station is obtained from the following 
equation (see eq. l) : 


where 
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When the indicated integration is performed, the moment relation for the 
constant-loading case becomes 
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For the first extreme -loading condition the loading varies linearly from 
the root chord to tip chord with the tip section loading 50 percent less 
than the root section loading. The total load is the same as for the 
constant -loading case. The relation of the loadings is expressed as 



r x 

p,c’ dy' = L = / pc’ dy' 
J 0 


where under the assumed loading condition 


(A5) 



(A 6) 


Substituting equation (A 6) into equation (A5) and then integrating gives 
the following relation 



Solution of this equation for the loading constant Kj yields 
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When the relation (A 6) is substituted into equation (A3) and then inte- 
grated, the moment is obtained for the first extreme- loading 

condition 
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The expressions for and % were substituted into equa- 

tion (A2) and the indicated integration performed graphically for taper 

(dOg) 

ratios of zero and unity. Curves showing the values of — along 

Ax b 

the span for these two taper ratios are presented in figure 17 . 

For the second extreme -loading condition the loading varies linearly 
and the r tip-section loading is 50 percent greater than the root-section 
loading. The relation of the loading is expressed in equation form as 



cLy 1 



pc ' dy' 


(A9) 


where under the assumed loading condition 


P 2 



(A10) 


As for the previous condition the loading constant was determined and 
is as follows: 


K 2 = 


6 + 6x 

7 + 8x p 


(All) 


When the relation (A10) is substituted into equation (A3) and then inte- 
grated, the bending moment Mgg is obtained for the second extreme con- 
dition of loading 
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The expressions for Mgg and Mg were substituted into equa- 
tion (A2) Pl ace 0 f and indicated integrations per- 

formed graphically for taper ratios of zero and unity. Curves showing 


the values of 


(**b)2 

Ax b 


along the span are also presented in figure 17 . 


I ‘-"Bn* l/-"BJp 

The variations of and are not presented outboard of 

the 0.8-semispan station. Near the tip the variations obtained by the 
method used in this appendix are not significant because the assumed 
loading and the variation of skin thickness in this region are not 
representative of the practical case. The conditions near the tip were 
modified for the constant-loading case used in the present report in 
order to account for these effects (-see discussion of angle-of -attack 
change due to bending in the section entitled "Method of Analysis"). 

If similar modifications were made for the cases of extreme loading 
considered here, the variations shown in figure 17 would approach a con- 
stant value (zero slope) at the. tip for finite taper ratios. The cal- 
culated variations are considered realistic inboard of the 0.8-semispan 
station. In this range the variations presented show that the effect 
of load distribution is greater for the case of zero taper ratio than 
for the case of unit taper ratio and, as would be expected, is greater 
near the tip than near the root. The magnitude of the effects of these 
alterations in load distribution from that assumed in the present paper 
(20- to 25-percent change in for outboard stations) are large 

enough to affect the preciseness of the quantitative results obtained by 
the method of analysis used in the present report but would not appear 
to affect the qualitative results. 
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Figure 1.- Typical, swept wing shcr 


used in the rireBent analysis. 
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(a) Plan-form taper ratio 0; thickness taper ratios 1.0 and 0.5 
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(b) Plan-form taper ratios 0.25 and 0.50; thickness taper 

ratios 1.0 and 0.50. 


Figure k.- Continued 
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(c) Plan-form taper ratios 0.75 and 1.0; thickness taper 

ratios 1.0 and 0 . 50 . 

Figure 4.- Concluded. 
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Figure 5*- Typical spanwise variations of the change in the angle of 
attack produced by wing bending for various plan-form and thickness 
taper ratios . 
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Figure 6.- Typical spanvise variations of the change in lii't produced 
by wing bending for various plan- form anfl thickness taper ratios. 
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Figure 10.- Variation with angle of sveepback of "bending stress per g normal acceleration, 
ratio of lift-curve slopes of flexible and rigid wing, and shift in aerodynamic center due 
to bending showing effects of aspect ratio, sweepback angle, and thickneBB ratio. Duralumin 
wings ; plan-form taper ratio 0.5; constant thickness ratio along span; symmetric parabolic- 

section; structural weight — W -- =* 0.2; Mach number, 0.9j altitude, 30,000 feet. 
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Figure 11.- Variation with aspect ratio of bending stress per g normal acceleration, ratio of 
lift-cufve slopes of flexible and rigid wing., and shift in aerodynamic center due to bending 
for. various sveepback angles showing effect of plan-form taper ratio. Duralumin wings; con- 
stant section thickness ratio of 0.03; symmetric parabolic-arc section; structural weight 
— r— = 0.2; Mach number, 0.9; altitude, 30,000 feet. 
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Figure 12 .- Variation with angle of sveepback of stress per g normal 
acceleration, ratio of lift-curve slopes of flexible and rigid wing 
and shift in aerodynamic center due to bending for thickness taper 
ratios of 1.0 and 0.5. Duralumin wings; aspect ratio 4; thickness 
ratio 0.08 (average of the root and tip section thickness ratios); 
plan-form taper ratio 0.5; symmetric parabolic -arc section; struc- 
tural weight - W - - = 0 . 2 ; Mach number, 0.9; altitude, 30,000 feet. 
























Figure 16.- Variations with, aspect ratio of the center of pressure of 
the loss in lift due to wing, tending and the tail area required to 
position the point for neutral static stability with respect to 
angle of attach at the center of pressure of the loss in lift when 
a rigid tail and fuselage are assumed. Wing taper ratio 0.5; con- 
stant thickness ratio. 





